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ABSTRACT 

Aims. We present and analyze e Aurigae data concerning the evolution of the H ff line on the occasion of the 2009 International 
observation campaign launched to cover the eclipse of this object. 

Methods. We visually inspect the dynamical spectrum constructed from the data and analyze the evolution with time of the EW 
(Equivalent Width) and of the radial velocity. 

Results. The spectroscopic data reveal many details which confirm the complexity of the e Auriga; system. The object is far from 
being understood. In particular, according to our measurements, the eclipse duration has been underestimated. A complete analysis of 
details revealed by our data would require much time and effort. Observers are encouraged to continue monitoring the H„ line out of 
eclipse in the hope that it will provide further important information. 

Key words, e Auriga; - epsilon Aurigae - 2009-201 1 eclipse - eclipsing binary - spectroscopy - H„ line - equivalent width - radial 
velocity - amateur spectroscopic survey. 



1. Introduction 

e Aurigae is one of the most intriguing eclipsing star sys- 
tems which has puzzled astronomers for nearly 200 years. 
The main eclipsing period is close to 27.1 years and the 
first spectroscopic surveys were undertaken during the 
1929 and 1956 eclipses. A large campaign was also or- 
ganized for 1982-1984. For a review of literature prior to 
the 2009-2011 eclipse, see lGuiNAN & DeWarf 20081 There 
are also numerous papers being prepared as a result of the 
2009-2011 eclipse. Despite the concenuated efforts, some 
aspects of e Aurigae remains a mystery 

e Aurigae is classified as an A8Iab star in an Algol type 
eclipsing binary system (Simbad) The prevailing model 
is of an F type star with a hot clumpy Hydrogen disc 
and an object of unknown nature which produces an 
eclipse phenomenon lasting almost 2 years every 27.1 
years. There may be a mid-eclipse brightening but solar 
proximity makes the photometry suspect at those times. 
Recently, it has been suggested that the eclipsing object 
is a 550 K dusty disk seen edge on, heated on the side 
facing the F star to 1100 K which may contain a B5V 
star (IHopkins et al. 201TT) that could contribute to emission 
wings surrounding H a line. Light curves feature 0.1 mag 
variations both inside outside eclipse. Variations have also 
been observed in the Equivalent Width (EW) and Radial 
velocity of spectral lines outside eclipse. These variations 
might be F star oscillations and wind. 



During 1982-84 eclipse, this star was studied by ama- 
teur observers using multiband photometric methods. The 
e Aurigae system was not clearly described despite all the 
acquired data. 

Twenty-seven years later, an international campaign 
was organized to manage both spectroscopic and photo- 
metric observations by amateur observers with the aim of 
producing data with improved time resolution compared 
with that achieved during previous eclipses. In this article, 
we present amateur spectroscopic H a line monitoring from 
February 12th 2008 to November 12th 2011. 

2. Observations 

In 2008, Jeff Hopkin^] organized the International observa- 
tion campaign of the 2009 e Aurigae eclipse. We acquired 
more than 250 high resolution spectra of the H a line cover- 
ing the three years around eclipse. These show significant 
variability throughout this period. The effect of the eclipse 
is clearly seen in this line from the end of April 2010 to end 
of April 2011. The spectra were recorded by five observers 
in Europe. 

Most observations were made using LHIRES30 spec- 
trographs. C. Buil used an eShell spectrograph that covers 



1 http : //www . hposof t . com/Campaign89 . html 

2 LHIRES3 and eShell are products from Shelyak Instrument, 
Grenoble, France: http : / /www . shelyak . com 
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wavelengths from 4500 A to 7000 A. Telescope diameters 
were between 0.2 m to 0.3 m. Spectral resolution is above 
10,000 and most of the time around 15,000. Mean expo- 
sure time was 2, 000 s. All setups are reported in table [T] 

3. Reduction and analysis method 

The raw observations are available from Robin 
Leadbeater's e Aurigas survey web paged. The spec- 
tra were reduced using standard procedures to produce 
calibrated and normalized lines profiles. Most of the 
reduction was done using SpcAudacd3 pipelines. The 
reduction steps were: geometric corrections, preprocessing 
and registration. Then line profiles were extracted with 
sky background subtraction. Wavelength calibration was 
done using calibration lamp spectra before and after each 
acquisition series to minimise the effects of calibration 
drifts. The instrumental response was then removed. An 
offset was computed using telluric lines to achieve a final a 
wavelength calibration RMS uncertainty of 0.03 A. Finally 
a heliocentric velocity correction was applied depending 
on the observation date. All wavelengths A are given in A 
(Angstrom). 

Equivalent Width measurements were computed be- 
tween /16550 and A6511 using linear integration and an ex- 
tracted continuum obtained from a fit to the local star con- 
tinuuirjfl The Chalabaev algorithm (IChalabaev 19831) was 
used to estimate the uncertainty, which is mostly depen- 
dent on the signal/noise ratio, and appears to overestimate 
the uncertainty compared with the actual scatter observed 
around the long term trend. 

Radial velocities were computed in two steps because 
the Hq, line is asymmetric: 

1. The Gaussian flank of the line was reproduced on the 
opposite side of the symmetry axis; 

2. A Gaussian fit of these two flanks gives a measure of 
the line center. 

A dynamical spectrum was computed using 177 spectra 
corrected to heliocentric velocity and cropped to ,16550- 
/16575. A linear interpolation was used to produce an 
image with a 1 day sampling interval. Such interpola- 
tion doesn't introduce bias for our analysis as the pur- 
pose of Fig. [6] is to show global behavior of the eclipse 
spread over several hundred days. Dates are logged in MJD 
(ie. JD-2400000). All computations were performed in 
SpcAudace. The monitoring covers a period of 719 days. 
Most of the information generated by our monitoring of the 
H a line is contained in this dynamical spectrum. 

3 http : //www . threehillsobservatory . co . uk/ 

4 http : //bmauclaire . free . f r/spcaudace/ 

s About EW's computation : http : //bmauclaire . free . fr/ 



Analyzing this complex image turns out to be cumber- 
some however. This is the reason why we have simplified 
the analysis by concentrating on the evolution with time of 
the EW that can be compared to V. mag, and of the radial 
velocity that can be used to study eclipsing phenomena. Of 
course we have to keep in mind that EW loses its physical 
meaning when applied to complex line profiles lines that, 
as it is the case for e Aurigae , are likely to result from a 
combination of several sources. But before analyzing these 
quantities, let us first examine the spectral line profiles at 
dates that show important transitions. 

4. Behaviour of the wings 

Outside eclipse, the H a line profile comprises a central ab- 
sorption core flanked by emission features on the red and 
blue wings (see Fig. [T]). These features are highly variable 
as described by Golovin (IGolov in 2008). 



epsAur- 12.881/12/2008 - SC12p + LHIRES3 2400 g/mm - 0.1 1 5 A/pixel - B. Mauclaire - 4x600 




6550 6555 6560 6565 

Wavelength (A) 



Fig. 1. H a wings out of eclipse phase. 

In the region of the H a line are absorption lines iden- 
tified as telluric lines at /16543.91, /16547.71, ,16548.32, 
,16552.63, -16557.17, /16568.81, ,16572.09, ,16574.85 and 
,16586.68. 

As we can see in Fig. |2] and |6j from MJD 55250 
onwards there was additional absorption in the core 
which broadened rapidly engulfing first the red emission 
and by MJD 55340 also the blue emission component. 
Note this is in contrast to the KI ,17699 line absorption 
which started decreasing in intensity during this phase 
(ILeadbeater et al. 201 II) . 

During ingress and into totality (see spectra at MJD 
55390.44 and MJD 55496.42), through the mid eclipse 
point, the absorption core deepened and broadened slightly 
on the red side. The additional absorption then moved to 
the blue and at MJD 55520 the red emission feature reap- 
peared. 

At the beginning of the decreasing phase (see spec- 
trum at MJD 55627.46), The H ff line became narrower with 
an emission component at the red side. Then from MJD 
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Table 1. Observation setups. 



Observer name 


Telescope 


Spectrograph 


Resolution 


Mean SNR 


Spectral range (A) 


Number of spectra 


C. Buil 


SCT 0.28 m 


eShell 


10000 


164 


4500-7000 


96 


T. Garrel 


SCT0.21 m 


Lhires3 2400 1/mm 


15 000 


108 


6500-6630 


87 


B. Mauclaire 


SCT 0.30 m 


Lhires3 2400 1/mm 


15 000 


224 


6520-6690 


34 


R. Leadbeater 


SCT 0.25 m 


Lhires3 2400 1/mm 


15 000 


98 


6500-6700 


25 


A. Lopez 


SCT 0.28 m 


Lhires3 2400 1/mm 


15 000 


164 


6500-6700 


5 



Time evolution of eps Aur from 2455250.403 to 2455878.39 
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Fig. 2. H a wings at key dates. 

55853.42, the blue edge emission component returned as 
just before eclipse. After the main eclipse phase, (see spec- 
trum at MJD 55878.39), the blue and red wings were both 
present but small, starting to resemble the before eclipse 
profile (Fig. CD- 

At the end of the survey period there still appears to be 
an excess absorption on the blue side of the central absorp- 
tion region compared with typical pre-eclipse spectra, pos- 
sibly due to the continued presence of the eclipsing disc. 
However the inherent variability of this at all phases makes 
the statement uncertain. 

Today's understanding (1Stencel2011I) is that the F star 
is semi-stable and capable of producing variability in lines 
in and out of eclipse. The disk is only modifying the optical 
spectrum during its passage. 

5. Equivalent width evolution with time 

EW measurements were computed between /16550 and 
A6511. Although the signal-to-noise ratio varies between 
observations and includes telluric lines which impact on 
EW, the effect most of the time is rather small. The data 
quality allows a reliable estimation of the EW. As men- 
tioned earlier, the single quantity EW is a gross simplifica- 
tion of the complex nature of the line. 



However, this quantity is the integral of the distribution 
of luminosity versus wavelength. It can thus be compared 
to similar integrals such as the V magnitude. As shown in 
Fig. and |U equivalent width (EW) and V magnitude (V 
mag.) are anticorrelated. Given that EW>0 for absorption 
lines, the eclipsing object occults the F star Hydrogen disk 
as the first minimum in EW and in V mag. evolution are 
both close to MJD 55250 and second minimum are also 
both close to MJD 55630. These dates define totality inner 
limits (see Table 2). 



Main contact 


MJD (day) 


first 


55070 + 3 


second 


55250 ± 2 


third 


55630 + 2 


fourth 


55800 ± 3 



Table 2. Contact dates during eclipse. 

While second and third contacts (and mid-eclipse) 
times are well-defined by our EW vs Date trend (Fig. [3]), 
the definition of first and fourth contact times are less ob- 
vious here and do not correspond to photometric contacts: 
these dates are likely to be linked with the densest ends of 
the eclipsing object. 

Equivalent width EW (6550-6577 A) variations within time 




55200 55300 



55500 55600 



Date (JD-2400000) 

Fig. 3. EW evolution within time. 

During the eclipse phase and outside it too, there are 
many small variations in EW and V mag. This suggests 
that the occulting object and F star Hydrogen disk may be 
clumpy. The F star may have also an intrinsic pulsating 
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Fig. 4. V mag. evolution within time. Curtesy AAVSO. 



activity (Kemp 19861 and IS tencel 20111) that produces such 
variations. 

Hq. EW has irregular variations like small steps during 
its increasing and its decreasing phases. Similar behaviour 
has been observed in the KI /17699 absorption line. It has 
been interpreted as an indication of structures (possibly 
ring like) within the disc (Leadb eater & Stencel 20101) . 
Continued observation during egress may help to clarify 
this. 

6. Radial velocity evolution with time 

Fig. |2] and [6] show how the line varied in radial veloc- 
ity. During eclipse (see spectra at MJD 55390.44 and 
MJD 55496.42), the absorption line became red shifted 
(+14.79 ± 1.37 km/s). During the end phase of the eclipse 
(see spectrum at MJD 55627.46), the H a line first returned 
to the position seen at MJD 55390.44 and then the absorp- 
tion line became blue shifted (-31.59 ± 1.54 km/s at MJD 
55853.42). See Table 3 for measurements at key dates. 



MJD (day) 


Radial velocity (km/s) 


55199.24 


+ 19.46 ± 1.54 


55390.44 


+ 14.79 ± 1.37 


55436.49 


-38.10 ± 1.54 


55521.46 


-60.36 ± 1.54 


55819.46 


-31.59 ± 1.54 



Table 3. Radial velocity measurements at key dates. 



An emission component (Fig. [5]) appeared in the core of 
the Hq. line close to the rest wavelength from MJD 55150 
onward as the absorption increased in this region. This be- 
came more clearly defined as the surrounding flux level 
dropped further and moved across the region from red to 
blue. The shape of the emission component is revealed 
as the absorption region broadened and swept across it 
through mid eclipse. It is clear that the constant emission 
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Fig. 5. Emission component at H a absorption line bottom. 

component is only revealed as the surrounding flux level 
drops. 

During our survey, the emission component measured 
by by gaussian fitting, remained centered on 6562.71 ± 
0.03 A, which, in terms of radial velocity, amounts to 
-4.97 ± 1.54 km/s. This does not account for the e Auriga? 
systemic radial velocity estimated to -2.26 ± 0.15 km/s 
(IStefanik 20101) . leading to a correction of +0.049 A. Note 
that Fig. [6] is not well enough resolved to see such a small 
shift. The computation was done on the line profiles. 

7. In quest of new models 

The light variations of Aur both in and out of eclipse have 
been the object of many studies. 

During the 1983 eclipse Kemp and Henson ( IKemp 1 986) 
analyzed polarization data. They suggested that the F-star 
is a non-radial pulsator and that its surrounding disk is 
tilted with respect to the orbit. 

In 1991, Ferluga suggested ( IFerluga 19911) that to ex- 
plain the shape of the light curve, the disk is not a contin- 
uous aggregate of dust, but instead a series of rings with 
a Cassini-like division. This model was more or less vali- 
dated by observations. 

During this eclipse, a wide variety of observations have 
been undertaken: infrared, ultraviolet, interferometry, pho- 
tometry, high resolution spectral monitoring. Thus a con- 
siderable amount of information (see IStencel 20101 and 
Hopkins et al. 2011] for an overview) is now available. It 
now remains to develop a model that fits all the data to 
hand. Undoubtedly the high resolution spectral monitoring 
data will be very important for constraining these models. 

8. Conclusions 

Amateur spectroscopists are now able to monitor bright 
targets with a spectral resolution of about 15, 000. Suitably 
equipped amateurs constitute a team with long term mon- 
itoring capacity which is widely distributed over the 
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Fig. 6. Dynamical spectrum of H„ line from MJD 55120.46 to 
MJD 55878.39. Interpolation between spectra was used to get a 
smooth image. Sinusoidal lines at both H a line edges are telluric 
lines lying in lines profiles that are corrected from heliocentric 
velocity. The black dashes along the right axis show the actual 
observation dates. 



troscopic contributions to the campaign over the coming 
months and years, especially those covering the Na D lines. 

We have observed similarities and discrepancies be- 
tween the EW and V magnitude evolution with time. The 
discrepancies remain to be explained which is beyond the 
scope of this letter. We also were able to define key dates 
in the eclipsing phenomenon. However much remains to 
be analyzed. Obviously the Hq, monitoring brings a lot of 
information which should place many constraints on the 
models conceived by scientists about e Aurigae. 

In any case we hope that this information will help sci- 
entists to solve the mysteries hidden behind this fascinating 
object. 

Acknowledgements 

This campaign would not have been possible without Jeff 
Hopkins's motivation and constant efforts to encourage 
amateur observations. The authors acknowledge Robert 
Stencel for his help and encouragements thoughout this 
campaign. 

All of the data presented in this letter were obtained 
from amateurs observers. The authors acknowledge also 
Eric Barbotin, Stephane Charbonnel, Valerie Desnoux, 
Stanley Gorodenski, Keith Graham, Torsten Hansen, 
James Edlin, Brian E. McCandless, Eric Sarrazin, Lothar 
Schanne, Jose Ribeiro, Francois Teyssier, Olivier Thizy 
and John Strachan for all the work they carried out on this 
mysterious star. 



Chalabaev, A. and Maillard, J.P. 1983, Astron. & Astrosphys., 127, 279. 
Ferluga, S. and Mangiacapra, D. 1991, Astron. & Astrophys., 243, 230. 
Golovin, A., Kuznyetsova, Y. and Andreev, M. 2007, Odessa Astron. Publ., 20, 
55. 

Guinan, E. F., Dewarf, L. E. 2002, in Exotic Stars as Challanges to Evolution, 

eds. C.A. Tout and W. Van Hamme, ASP Conf. Ser. 279 Astron. Soc. 

Pacific, San Francisco, 121. 
Hopkins, J., Stencel, R., Leadbeater, R., 2011, Epsilon Aurigae Eclipse 

International Campaign Newsletter, No. 24 (Fall/Winter). 
Kemp, J., Henson, G., Kraus, D., Beardsley, I., Carroll, L., Ake, T., Simon, T. 

and Collins, G. 1986, Astrophys. J., Lett. Ed., 300, LI 1. 
Leadbeater, R„ and Stencel, R. 2010, http://arxiv.org/pdf/1003.3617v2.pdf (18 

March). 

Leadbeater, R., Buil, C, Garrell, T., Gorodenski, S., Hopkins, J., Mauclaire, 
B., Ribeiro, J., Schanne, L., Thizy, O., and Stencel, R. 2011, Bull. Amer. 
Astron. Soc, 43, 2 57.04. 

Simbad database. 201 1, e Aurigae web page 

Stefanik, R. P., Toress, G., Lovegrove, J., PerAa, V. E., Latham, D. W., Zajac, J., 

and Mazeh T. 2010, Astron. J., 139, 1254. 
Stencel, R. 2010, "Epsilon Aurigae in Total Eclipse 2010, Mid-eclipse report", 

Univ. Denver Obs. (http://arxiv.org/pdtyi005.3738.pdf) 
Stencel, R. 2011, private communication (December 7). 



References 



planet. Our H a monitoring of e Aurigae shows that, con- 
trary to what was forecast, the effects of the eclipse ex- 
tended beyond December 2011. Post-eclipse observations 
are needed. R. Stencel welcome any outside eclipse spec- 
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